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ABSTRACT 


A method for calculating the average spatial 
and energy characteristics of hadron-lepton 
cascades in the atmosphere is described. The 
results of calculations for various strong 
interaction models of primary protons and 
nuclei are presented. The sensitivity of the 
experimentally observed EAS characteristics 
to variations of the elementary act parameters 
is analyzed. 


The theoretical analysis of hadron-lepton cascade 
propagation in the atmosphere is possible at present by 
two methods: a) Monte-Carlo cascade simulation; b) solving 
the set of kinetic equations. Despite the advantages, the 
Monte-Carlo method is at the same time the most cumbersome 
and demands much computer time. The main drawback of the 
second method is that shower characteristics for fixed 
primary energy are obtained, whereas in the experiment the 
showers with fixed Ne are studied. However the comparison 
of the results obtained by both methods shows that the de- 
pendence of the average shower characteristics on the 
strong interaction parameters turns out to be qualitative- 
ly the same, while the quantitative difference is small. 

As an example of this, we may point out the conclusion that 
the average experimental characteristics of EAS can be de- 
scribed using the strong interaction model with increasing 
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cross sections, growing plateau in the pionization region 
and asymptotic scaling in the fragmentation region; this 
conclusion, being made from the solution of diffusion 
equations /I/, was confirmed by the Monte-Carlo simulation 


results / 2/. 

In the present paper we describe a method of solving 
the three-dimensional set of kinetic equations for shower 
components, based on the calculation of the" contributions 
from the successive generations by Monte-Carlo evaluation 
of the corresponding multidimensional integrals. 

The integral spectrum of the (i + 1 )-th pioh generation 
in EAS for depth Z and primary energy E 0 can be present- 
ed in the form: 






dV 


i+ < 


Ei + i 
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in*. 



where 7T(>E, £) , TT(p,2.) and N(p,£) are the inte- 

gral and differential momentum spectra of plons and nuc- 
leons at depth Z, and (p, ft ) and p 0 ) 

are the invariant inclusive spectra of pion production in 
NA and TT A interactions. The expressions for the nucleon 
and muon spectra can be written in the same manner. For 
Monte-Carlo computation of integral (1) a set of interac- 
tion depths in the atmosphere and particle momenta are ran- 
domly simulated. From the known momenta and depths, we cal- 
culate the diatance from the shower axis to the particle's 
arrival point at the observation level. The value of the 
integrand is calculated for the simulated. set. Further on, 
knowing the energy and the distance of the particle from 
the axis, we can obtain the integral and differential spec- 
tra over E and r . The spatial distribution function 
(SDF) of electrons is obtained multiplying the ft -produc- 
tion spectrum r(E,H, r) by the SDF for photon initiated 
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ahowera in the approximation of ref./3/. The accuracy of the 
calculation of the integral depends on the variance and the 
number of simulations of the integrand values. The integral 
was smoothed to minimize the dispersion. The method will be 
described in detail in a forthcoming paper. 

This work presents the results of calculations for the 
models of multiple production. The first one is the so-call- 
ed "standard' 1 model of the scaling type, where the inclusive 
spectra are taken in the Hillas parametrization /4/, and 
^ a 0.37 GeV. The second is the model of quark— gluon jets 
/5/ , where the target nuclei effects are taken into account 
according to ref./6/. In this model, as is well known, there 
takes place the cross section rise, the plateau grows and 
the scaling is weakly violated in the fragmentation region. 
In this model the spectra differ quantitatively from those 
used earlier, mainly in the fragmentation region. The mean 
transverse momentum in this model increases logarithmically 
with energy. 



Figure 1 shows the radial 
distribution of -T 2 "/ N e , 
where is the SDF at 

p 

700 g/cm . The figures near 
the curves denote the pri- 
mary energy in GeV. S labels 
the "standard" model curve, 

Q - the models of the quark- 
gluon jets. The curves dif- 
fer insignificantly, since 
the SDF of EAS is mostly due 
to the electromagnetic in- 
teractions. However in the 
model a more rapid energy 
dissipation takes place and 


the SDF is a bit wider. 

Figure 2 presents hadron and muon distribution in a 
circle of radius with E > 5 GeV. The same is shown in 
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Pig. 3 for an energy cut-off 1 TeV. The small energy par- 
ticles number and the SDP show weak moel dependence. As ex- 
pected, the high-energy hadrons are concentrated in the 
narrow region near the shower axis. The same is true for 
high-energy muons. The SDP of high-energy hadrons is the 
most model-dependent characteristic. In the Q model the 
total number of muons with energy > 5 GeV is greater than 



in the 5 model , whereas the number of muons with E > 1 TeV 
is almost the same but the SDP in the Q model is wider which 
is due to both rise of the cross section and the average 
transverse momentum. The total number of hadrons with E > 1 
TeV is larger in the S model as compared to that in the Q 
model by a factor of 2.5 in the range of E c and T . 

The authors are indebted to the participants of the 
ANI-84 workshop for useful discussions. 
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